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Observation of ferroelectricity among non-d° systems, which was believed for a long time an unrealistic 
concept, led to various proposals for the mechanisms to explain the same (i.e. magnetically induced 
ferroelectricity) during last decade. Here, we provide support for ferroelectricity of a displacive-type possibly 
involving magnetic ions due to short-range magnetic correlations within a spin-chain, through the 
demonstration of magnetoelectric coupling in a Haldane spin-chain compound Er 2 BaNi0 5 well above its 
Neel temperature of (T N =) 32 K. There is a distinct evidence for electric polarization setting in near 60 K 
around which there is an evidence for short-range magnetic correlations from other experimental methods. 
Raman studies also establish a softening of phonon modes in the same temperature ( T) range and 
T-dependent x-ray diffraction (XRD) patterns also reveal lattice parameters anomalies. Density-functional 
theory based calculations establish a displacive component (similar to d°-ness) as the root-cause of 
ferroelectricity from (magnetic) Ni0 6 chain, thereby offering a new route to search for similar materials 
near room temperature to enable applications. 
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The field of magnetoelectric (ME) coupling and multiferroicity involving two correlated phenomena, namely, 
ferroelectricity and magnetic ordering, bearing relevance to applications, has expanded recently. However, 
this kind of multiferroicity has been considered rare in transition-metal systems, as magnetism involves 
partially filled d orbital and ferroelectricity requires empty d orbital. The observation of magnetism-induced 
ferroelectricity in some materials led to the proposals of several mechanisms 18 , for example, lone-pair stereo- 
chemistry, role of spin-spiral structure, charge ordering, and geometrical considerations. However, a displacive- 
type mechanism involving off-centering of the magnetic ion proposed in Ba-doped SrMnCV' 10 is uncommon 
among stoichiometric (undoped) compounds. Here, we demonstrate that the one-dimensional short-range 
magnetic correlation plays an important role to induce ferroelectricity above long-range magnetic ordering in 
an undoped spin-chain compound possibly by this mechanism. 

The compound under discussion is a Haldane spin-chain compound, Er 2 BaNi0 5 , crystallizing in an ortho- 
rhombic structure (Immm). We show that this family, R 2 BaNi0 5 (R = Rare-earth), is characterized by a strong 
ME coupling and potentially prone for multiferroicity due to a displacive-type mechanism. Though this family of 
insulating oxides aroused considerable interest due to Haldane spin-chain magnetic anomalies 1122 , very little 
work has been reported 17 ' 21,22 to understand dielectric or multiferroic behavior. Experimental studies on a few 
members of this family e.g. Ho 17 and Dy 22 were reported which revealed ferroelectric behaviour, but the micro- 
scopic mechanism of the observed behaviour has not been clear. Keeping this in mind, we have explored ME 
coupling on the Er member through macroscopic and microscopic measurements, augmented by Raman scatter- 
ing and XRD studies, and electronic structure calculations. The orientation of the magnetic moments has been 
known to be decided by the anisotropy of the rare-earth ions and therefore the primary motivation with which 
this study was undertaken was to explore whether this factor anyway controls appearance of multiferroicity in this 
series. For instance, in Dy 2 BaNi0 5 , in which long-range antiferromagnetic order sets in at (T N =) 58 K, the 
magnetic moments of Dy and Ni at 1.8 K lie almost along c-axis, with these moment- directions getting gradually 
decoupled with increasing temperature rotating towards a-axis (along Ni-O-Ni chain axis). On the other hand, in 
contrast to other R members, in the Er case, the directions of Er and Ni moments point along a-axis at all 
temperatures below T N 13 , thereby suiting to the motivation of this work. 
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Figure 1 | Dielectric constant and the loss factor as a function of 
temperature measured with a frequency of 50 kHz in the presence of 
external magnetic fields are plotted in (a) and (b) for Er 2 BaNi0 5 . The 

plots of tand shown in (b) are to be used primarily to highlight the feature 
near 45 K in the presence of H. The heat-capacity data ( C/T vs T) measured 
under different magnetic fields are plotted in (c); in inset of (c), the data is 
plotted in a wider temperature range to highlight the weak feature around 
T N (marked by an arrow) and a broad anomaly above T N . 

Results 

The T-dependence of dielectric constant (e') and loss factor (tand) in 
the presence of various magnetic fields (H) are shown in figure 1 (a,b) 
for a frequency (v) of 50 kHz. The behavior of heat-capacity (C) is 
also shown in figure 1, primarily to understand the features in com- 
plex dielectric data. In the C/T curve, in the absence of a magnetic 
field, there is a weak jump at 32 K at the onset of magnetic ordering 
(see inset of Fig. lc), followed by a gradual fall with decreasing 
temperature till about 15 K, below which there is another feature 
manifesting as an upturn and a peak at 8 K. It is clear from the 
figure lc that the temperatures at which the upturn and the peak 
occur (below 20 K) get lowered gradually with increasing magnetic- 
field (H). A broad hump is also observed spreading over a wide T- 
range around 40 to 60 K, which could be attributable to (intrachain) 
short-range magnetic order, also considering that the presence of 
magnetic short-range order above T N was reported by a spectro- 
scopic study in this compound 23 . Keeping these observations in 
mind, we turn to dielectric data. The values of tand are very low 
consistent with the insulating behavior of this compound. We noted 
a fall in e' below 8 K for 5 kHz with a concomitant anomaly (a peak) 
in tanS and this feature was found to exhibit a v-dependence with 
this fall appearing at —10 K for 70 kHz (not shown here). Similar V- 
dependence was observed in ac susceptibility as well (not shown 
here) and we will not focus this aspect as this is not the aim of the 
article. Though it is not easy to infer a change at T N from the plots of 
£.'(T) (see, for instance, the curve shown in Fig. la for v = 50 kHz) 
obtained in the absence of a magnetic field, a careful comparison with 
the plots taken in the presence of H reveals an anomaly around this 
temperature. In e'(T}, the curves for different H values essentially 
overlap above —40 K, and these tend to bifurcate clearly below a 
temperature close to 32 K, thereby offering a support for ME coup- 
ling. The corresponding feature observed in figure lb (for tan<5 (T)) is 



more intriguing, as it reveals a distinct separation of the curves mea- 
sured under different H at a temperature much higher than T N , that 
is, below —45 K, though the observed change in the value is very 
small near 45 K. We attribute this to a magnetic precursor effect 
before the long range magnetic ordering sets in, that is, to the exist- 
ence of short-range correlations proposed above. The value of e' 
tends to increase with H below T N . It is worth noting that the tem- 
perature at which a sudden drop (for the lOK-feature in figure la) 
occurs decreases with increasing H, for instance from 10 K for H = 0 
to about 8 and 4 K for 50 and 140 kOe respectively for 50 kHz, with 
this temperature closely tracking the trend observed for the corres- 
ponding feature in the C(T) curves (Fig. lc). The decrease of anom- 
aly/peak position towards low temperature with increasing magnetic 
field indicates dominant antiferromagnetic component around this 
temperature. The observed correlation near 10 K between dielectric 
and heat capacity results provides an additional support for the 
proposal of ME coupling in this compound. This dielectric behavior 
with the application of magnetic field is the same for other frequen- 
cies (5-100 kHz) as well, though the absolute values are different (not 
shown here). 

We have also obtained further evidence for ME coupling effect on 
the basis of isothermal magnetization (M) behavior and e' as a func- 
tion of H at different temperatures, the rate of variation of H being 
70 Oe/s (see figure 2). At 2 K, we see a distinct evidence for two 
jumps in M(H), a stronger one near 20 kOe and the other weak 
one near 30 kOe (Fig. 2a). Such metamagnetic transitions were 
reported at 4 K as well by Chepurko et al 23 . To bring out ME coup- 
ling, we show Ae' as a function of H at 2 K in figure 2b, where Ae' = 
(e'h = o- e'h)/ s'h = o- We find that there is a sudden upturn in the 
value around 20 kOe with another weak shoulder around 30 kOe, 
followed by a gradual non-linear rise at higher fields, similar to M(H). 
This correlation between M(H) and Ae'(H) data establishes ME 
coupling in this compound. 

The temperature dependent remnant polarization (P) was mea- 
sured during warming in the absence of electric-field after cooling the 
sample in an application of electric-field (E) of 220 kV/m under zero 
magnetic field and also in the presence of different H for the field- 
cooled condition (see the method section for details). The results 
obtained are shown in figure 3a. Interestingly enough, though long 
range magnetic order sets in at 32 K, remnant polarization sets in at a 
much higher temperature (—60 K) compared to T N . The sign of 
polarization was also found to get reversed with a reversal in the sign 
of electric field (inset of figure 3a), which is sufficient to claim the 
ferroelectric behaviour of the sample and similar ferroelectric beha- 
viour has been reported for Ho 2 BaNi0 5 17 and other materials also 1,2 . 
We endorse a possible explanation in terms of short-range low- 
dimensional magnetic interactions well above T N 16 , which in turn 
facilitates ferroelectric behaviour through magnetoelectric coupling. 
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Figure 2 | Isothermal magnetization (M(H)) and change (Ae') in 
dielectric constant, As' = (e' h = 0 - e' h )/ e'h = o f° r Er 2 BaNiO s . 
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Figure 3 | For Er 2 BaNi0 5 , (a) electric polarization, P, as a function of 
temperature (in the absence and in the presence of external magnetic 
fields), obtained as described in the text; inset: Pfor a change in the sign of 
electric-field for H = 0. (b) Pyroelectric current data; inset shows electrical 
resistivity data. 

It should be stressed that such reports bringing out multiferroic 
behavior due to isolated spin-chains are generally rare in the literat- 
ure and we have seen such effects in another spin-chain compound 
Ca 3 Co 2 0 6 recently 24 . As an evidence for ME coupling, the values of P 
are found to change with the application of magnetic-fields well 
below 60 K, in addition to a very small change in the ferroelectric 
transition temperature (around 60 K). The plot (see Figure 3b) of 
pyroelectric current, / pyro , also exhibits a peak around 55 K with its 
intensity and the peak-temperature varying with H marginally. It is 
worth noting that the overall variation of P in the magnetically 
ordered state is much larger than that observed for the Dy analogue 22 ; 
we observed a value of about 4 uC/m 2 for the Dy case by an applica- 
tion of 400 kV/m, whereas for the Er member, it is —100 (j,C/m 2 for 
an application of 220 kV/m only; for the Ho case, the results were not 
reported 17 for the absence of a magnetic-field to enable us to com- 
pare. It is to be noted that the P value almost saturates below mag- 
netic ordering T N . 

It is worth stating that the dc resistivity of this compound at room 
temperature is of the order of giga-ohm-cm and increases exponen- 
tially with decreasing temperature; below 120 K, it increases drastic- 
ally further (inset of figure 3b) and, below 90 K, it was not possible to 
measure with proper resolution due to the limitations of the electro- 
meter. This clearly establishes that this compound is a highly insu- 
lating system and whatever ME coupling we report here is not due to 
leakage current. 

We have also measured the complex impedance as a function of 
frequency at different temperatures applying different bias voltages 
(0.1-2 V) using the same LCR meter and the experimental set up 
which was used for dielectric measurements. It was noted that the 

current-voltage phase angle is 90° (above —89.5° below around 

70 K within the measured frequency range) for a wide frequency 
range (between 2-500 kHz). This means that the compound is 
undoubtedly a very good capacitor in this range, which is a pre- 
requisite for the investigations for the aim of this article. There is 
no change of feature or value with different bias voltages below 
100 K. Therefore, we are confident that our conclusions from dielec- 



tric measurements are intrinsic, without any dominant influence of 
extrinsic effects (like contact problem, grain boundary etc.). 

We have also obtained spectroscopic evidence for mode-softening 
attributable to the formation of polar state through Raman scattering 
experiments. Er 2 BaNi0 5 forms in the space-group Immm at room 
temperature 12 . The compound has 9 atoms per unit cell and, out of 
the 24 vibrational modes, only the Er and the O atoms participate in 
the Raman modes, while all others participate in the infra-active (IR) 
modes only 25 . Figure 4a shows vibrational spectra for representative 
temperatures. Certain modes marked by shaded bands (B lt B 2 , B 3 ) 
were found to lose intensity at lower temperatures, while others 
persist. The persisting modes that are marked as S 3 , S 7 , S 8 are 
assigned the symmetry of B lg> B 3g , and A g respectively and are 
Raman active. The modes S 5 and S 7 appear to be combinational 
Raman scattering modes. The mode assignment was done by com- 
paring the present spectra with those for other compounds of the 
Haldane chain family 25-27 . We assign the spectral features in the 
shaded bands to the IR modes which arise from any kind of disorder 
in the Ni-O chain. The key observation is that there is a dramatic loss 
of intensity of these disorder-based modes with decreasing temper- 
ature. We have to necessarily attribute the loss of intensity of the IR 
modes to the gradual development of magnetic order within this 
chain (as inferred from the C(T) as well). The intensity of these 
magnetic-disorder induced modes is thus proportional to the mech- 
anical susceptibility of the Ni-0 chain. If the susceptibility is higher, 
it is more probable for thermal fluctuations to generate disorder 28 . 
The inverse of the intensity is thus related to the stiffness of the Ni-0 
chain. Figure 4b shows the inverse of the intensity of the disorder- 
induced modes as a function of temperature. It is apparent that there 
is a distinct inverted cusp around 45 K, as though there is an order- 
disorder phase transition of second-order type. The signature of this 
transition can also be seen in the variation of the peak position (oj,) of 
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Figure 4 | (a) Vibrational spectra at representative temperatures for 
Er 2 BaNi0 5 . The features corresponding to vibrational modes are shown 
with shaded bands. The disorder-induced infrared modes are within the 
shaded-band. The modes S 3 , S 6 , S 7 , S 8 and S 9 are Raman modes of single 
scattering events, while S 5 and S 7 are Raman modes that involve multiple 
modes; (b) the inverse of the integrated intensity of the modes within the 
shaded-band as a function of temperature. The dash-dotted line is a guide 
to the eyes, (c) Variation of peak position of the representative modes 
Si-S 9 as a function of temperature. The shaded region marks the 
temperature range from 60 K to 45 K where there is a distinct change in 
slope. 
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some representative peaks Si...S 9 (marked in Fig. 4a) as a function of 
temperature (see Fig. 4c). The value of daVoTclearly shows a change 
at 60 K and in the temperature range 45 < T <60 K (shaded region), 
it is positive. We associate this mode stiffening to the distortion of 
Ni0 6 octahedra. Below 40 K, the sign of the slope is negative. This 
continuous softening of the modes at low temperatures is related to 
the renormalization of the phonon frequencies due to the presence of 
strong spin-phonon coupling in this system. 

The results presented above provide distinct evidence for soften- 
ing of the modes in the T-region where short-range magnetic-order- 
ing-induced ferroelectricity tends to set in. The observation of 
ferroelectricity in this family is indeed puzzling, especially since the 
experimentally determined crystal structure as shown in Fig. 5(a) is 
centro-symmetric. The Ba and rare earth ions occupy high symmetry 
positions in the tetragonal unit cell. In order to throw more light on 
this aspect, electronic structure calculations as well as phonon calcu- 
lations were performed (for calculation details see in method sec- 
tion). The generalised gradient approximation 29 was used for the 
exchange-correlation functional, though we also performed calcula- 
tions including an onsite Hubbard U-like term on Ni to treat elec- 
tron-electron interactions in a mean-field manner. Various magnetic 
configurations were considered; however, the lowest energy was 
found for the one with antiferromagnetic one-dimensional network 
of Ni-0 chains running along the a-direction. The coupling between 
chains was found to be extremely weak and —0.4 meV. As a starting 
point, the phonon dispersions were calculated for the reported struc- 
ture 13 along various symmetry directions using the small displace- 
ment method. Soft phonon modes were explored and displacements 
were made along the eigenvectors corresponding to the soft modes 
and the internal positions were again optimised. In every case, we 
worked with the experimental volume, optimising only internal 



coordinates. A soft mode was identified to develop at T point 
(Fig. 5(b)). Using the displacements suggested by the soft mode, 
we arrive at a lower energy structure which is non-centrosymmetric 
and corresponds to the space group lmm.2. Each Ni atom is sur- 
rounded by six oxygens, though the connectivity of these Ni0 6 
motifs is present only in the a-direction through corner-shared oxy- 
gens. The primary displacement in the ferroelectric structure 
involves the Ni atoms moving towards a pair of oxygens in the bc- 
plane. As a result, a pair of Ni-O bond-lengths become shorter while 
the other pair are longer. These become 2.177 and 2.187 A (Fig. 5(c)) 
when the experimental antiferromagnetic structure is imposed, while 
these are equal to 2.142 and 2.228 A (Fig. 5(d)) when a ferromagnetic 
structure is considered. The Ni-O bond-lengths along the chain dir- 
ection, remain unchanged. There are small dipole moments assoc- 
iated with the Ba-O and Er-0 network also. The polarization is 
computed to be 0.39 uC/cm 2 for the ferromagnetic solution, while 
it is scaled down by almost an order of magnitude for the experi- 
mentally observed antiferromagnetic structure. Further, increasing 
the inter-chain separation artificially in our calculations, thereby 
making the chains more one-dimensional, we find an increase in 
the observed polarization. It is to be noted that polarization is absent 
if we do not introduce the one-dimensional magnetic interaction in 
our calculation. 

In order to understand the mechanism for the observed polariza- 
tion, we have examined the density of states. The system is found to 
be an insulator. Broadly, for any transition-metal atom in an octa- 
hedral environment, there is a lifting of degeneracy of the d-orbitals 
with the orbitals with t 2g symmetry at lower energies compared to the 
orbitals with e g symmetry. The reduced symmetry of the Ni0 6 octa- 
hedron in the present case (D 2 h), lifts the degeneracy further; how- 
ever, we can still identify the d-orbitals which have t 2g symmetry and 
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Figure 5 | (a) The projection of crystal structure of Er 2 BaNi0 5 along a-direction; (b) the phonon dispersions are shown for the symmetry 
directions T-X and T-L; and (c) and (d) show the Ni-O bond lengths in the ferroelectric states for antiferromagnetic and ferromagnetic structures 
respectively, as obtained in the calculations. 
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those which have e g symmetry and so we continue the discussion in 
terms of the t 2g and e g states. The basic energy level diagram can be 
understood as consisting of completely filled t 2g orbitals on Ni, and 
half-filled e g orbitals. As the bandwidth is smaller than the exchange 
splitting, the system is insulating. 

The calculations described above reveal that there are changes in 
the Ni-O bond distances. In order to gain knowledge on the influence 
of ferroelectricity on the lattice parameters, we have carried out x-ray 
diffraction studies as a function of T down to 10 K. In the powder 
diffraction patterns, we have not observed any additional reflections 
in the whole temperature range, other than those known for this 
compound. Rietveld refinement (see supplementary information) 
was carried out to capture any signature of the phase transition. 
Refined unit-cell parameters and volume are plotted in figure 6. It 
is evident from this figure that the cell-parameter 'b' indeed shows a 
clear upward deviation from the thermal contraction behavior with 
decreasing T in the vicinity of ferroelectric transition temperature 
(around 70 K), apart from the one close to T N . This deviation is well- 
outside the error bar [—0.0001 A]. The unit-cell parameters 'a and V 
also reveal a weak change of slope in the same temperature regions. 
Naturally, the unit-cell volume also shows a distinct anomaly. We 
have also explored whether there is any change in the lattice sym- 
metry following ferroelectric transition by Rietveld fitting of XRD 
patterns for two ferroelectric subgroups, Imm2 and 1222, of (room 
temperature) Immm space group. For Imm2, oxygen at 8 1 site [of 
Immm] splits into two 4 d [of Imm2] Wyckoff sites, resulting in 
increase in refinable position parameters. For 1222, all three position 
coordinates [8 k site] of oxygen atom can be varied during refine- 
ment. We found that the pattern is equally well-fitted with all the 
three space groups (see supplementary information). Due to possible 
very small changes in position coordinates across the phase trans- 
ition, the present Rietveld refinements can not conclusively deter- 
mine the true symmetry of ferroelectric phase. Such difficulties are 
not unusual when the lattice distortions are weak, as remarked in 
Ref. 30 (also see, Ref. 7 in Ref. 30). In any case, the present XRD 
results provide evidence for a displacive-type ferroelectricity as indi- 
cated by changes in cell-parameters across the transition. 
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for Er 2 BaNi0 5 as a function of temperature determined from x-ray 
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region around which lattice constants reveal some anomalies. 



Discussion 

It is important to note that the distortion observed in our calculations 
is similar to what has been seen in d° ferroelectrics, but rarely seen in 
systems with a finite d-electron count. From the details provided 
above, it is clear that this is not the same as the exchange-striction 
driven mechanism observed for Ca 3 CoMn0 6 31 . The natural question 
at this point is why the displacive ferroelectricity proposed here is 
rare. In the past, systems exhibiting non- symmetric degeneracy lift- 
ing first-order Jahn-Teller distortions 32 34 as well as band-insulators 
exhibiting second-order Jahn-Teller distortions 7,8-35 have been 
believed to be possible candidates for displacive ferroelectricity 
among finite d-electron systems. The electronic structure discussed 
earlier puts the system into a class of band-insulators. However, every 
band insulator system does not turn ferroelectric and the question 
arises why is that the present family is so susceptible to ferroelectri- 
city. In the case of the classic example of BaTi0 3 ferroelectric 36 , it is 
known that the significant energy-gain for the ferroelectric structure 
comes from the band-energy gain arising from shorter Ti-O bonds 37 . 
For Ni 2+ oxides, one finds from an earlier work 38 by some of us that 
the charge-transfer energy decreased as one moves from three- 
dimensional NiO to a two-dimensional network of Ni0 6 motifs 
and then to a one-dimensional one. A smaller charge transfer energy 
implies a larger band energy gain, and hence ferroelectricity seems to 
be likely in these one-dimensional oxides. 

In short, we establish that the Haldane spin-chain compound, 
Er 2 BaNi0 5 , which has been previously shown to undergo long-range 
antiferromagnetic order below 32 K, is characterized by magneto- 
electric coupling with ferrolectricity setting in above T N . The key 
conclusion is that this family of magnetic materials is prone for 
multiferriocity under favorable circumstances due to the 
'Displacive-type mechanism', that distorts the 0 6 octahedra, thereby 
lifting the point of inversion symmetry. This is not so-commonly 
encountered in undoped materials with non-3d° electronic config- 
uration. Interestingly, short-range magnetic correlations appear to 
be adequate to trigger ferroelectricity due to this mechanism. We 
hope this work triggers further work to search for similar materials 
near room temperature. 

Methods 

Polycrystalline form of the compound, Er 2 BaNi0 5 , was prepared as described earl- 
ier 16 starting with stoichiometric amounts of high-purity (>99.95%) oxides, BaC0 3 , 
Er 2 0 3 and NiO. X-ray diffraction study established single phase nature of the 
materials. 

The dc x measurements were carried out using commercial Superconducting 
Quantum Interference Device (SQUID, Quantum Design, USA) and isothermal 
magnetization behavior was also performed with the help of a commercial Vibrating 
Sample Magnetometer (VSM, Quantum Design, USA) at various temperatures. The 
features in M(T) and M(H) curve are in broad agreement with those reported in the 
literature 13,14 ' 6 , thereby ensuring the quality of the sample. 

Heat-capacity measurements were done by a relaxation method using commercial 
Physical Properties Measurement System (PPMS, Quantum Design). 

Complex dielectric permittivity measurement were carried out using Agilent 
E4980A LCR meter with a homemade sample holder coupled to the PPMS. 
Temperature dependent complex dielectric permittivity was measured for various 
frequencies (5-100 kHz) at 1 V ac bias during warming (0.5 K/min). Isothermal 
magnetic-field dependence of dielectric behavior was measured at different tem- 
peratures (below and above T N ) at different frequencies (10-100 kHz). The rate of 
change of magnetic-field in all these measurements (magnetization and dielectric 
both) is 70 Oe/sec. Remnant polarization was measured with Keithley 6517A elec- 
trometer (using the same set up like dielectric measurement integrated to PPMS) in 
Coulombic mode (with automated integration of I pyro ) as a function of temperature. 
An electric-field of 220 kV/m was applied at 80 K (above ferroelectric temperature) 
to align the electric dipoles and then the sample was cooled to 6 K; subsequently, the 
electric-field was switched off and the capacitor (sample) was shorted for sufficient 
time to completely remove the stray charges (if any). P(T) curve was then obtained 
during warming (2 K/min). We have also obtained these curves under different 
magnetic fields with the application of both H and E at 80 K and then the specimen 
was cooled down to the lowest measured temperature in the presence of both the 
fields; then polarization was measured in the same manner as discussed above, only 
keeping the magnetic field switched on till the end of the measurements. 

The dc resistivity was also measured in the two-probe method using Keithley 
6517A electrometer. 
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Raman scattering measurements were performed in the backscattering geometry 
using a triple grating Raman spectrometer (T64000: Jobin Yvon) equipped with a 
liquid nitrogen cooled charge-coupled device. The excitation source was the 
647.1 nm line of a mixed gas laser (Stabilite 2018: Spectra Physics). Low temperature 
Raman experiments were carried out using a continuous flow helium cryostat 
(Microstat He: Oxford Instruments). 

Electronic structure calculations were performed within a plane- wave imple- 
mentation of density functional theory within the Vienna Ab initio Simulation 
Package (VASP) code 39,40 . A k-points mesh of 4 X 4 X 4 were used and an energy 
cutoff of 500 eV was used for determining the maximum kinetic energy of the plane 
waves included in the basis. The polarization of the ferroelectric structure was 
computed using the Berry phase method as implemented in VASP. 

The x-ray diffraction patterns were taken at different temperatures in the range 10- 
300 K in a symmetric Bragg- Brentano geometry, with horizontally scanning 20-arm 
using Rigaku-make rotating anode x-ray generator operated at 12 kW. Lowering of 
temperature was achieved using a liquid-helium cryostat (Oxford Instruments, UK), 
which is directly mounted on a goniometer from Huber, Germany. The high quality 
of the data can be judged from the low and flat background with high peak to back 
ground ratio (which is more than 500). Rietveld refinement has been carried out using 
FullProf Suite package 41 . 



1. Khomskii, D. Classifying multiferroics: Mechanisms and effects. Physics 2, 20 
(2009). 

2. Rao, C. N. R., Sundaresan, A. & Saha, R. Multiferroic and magnetoelectric oxides: 
The emerging scenario. /. Phys. Chem. Lett. 3, 2237-2246 (2012). 

3. Seshadri, R. & Hill, N. A. Visualizing the Role of Bi 6s "Lone Pairs" in the Off- 
Center Distortion in Ferromagnetic BiMn0 3 . Chem. Mater. 13, 2892-2899 
(2001). 

4. Newham, R. E., Kramer, J. J., Schultze, W. A. & Cross, L. E. 
Magnetoferroelectricity in Cr 2 Be0 4 . /. Appl. Phys. 49, 6088-6091 (1978). 

5. Ikeda, N. et at Ferroelectricity from iron valence ordering in the charge -frustrated 
system LuFe 2 0 4 . Nature (London) 436, 1136-1138 (2005). 

6. Aken, B. B. V., Palstra, T. T. M., Filippetti, A. & Spaldin, N. A. The origin of 
ferroelectricity in magnetoelectric YMn0 3 . Nature Mater. 3, 164-170 (2004). 

7. Rondinelli, J. M., Eidelson, A. S. & Spaldin, N. A. Non-d° Mn-driven 
ferroelectricity in antiferromagnetic BaMn0 3 . Phys. Rev. B 79, 205119 (2009). 

8. Lee, J. H. & Rabe, K. M. Epitaxial-strain-induced multiferroicity in SrMn03 from 
first principles. Phys. Rev. Lett. 104, 207204 (2010) 

9. Sakai, H. et at Displacement-type ferroelectricity with off-center magnetic ions in 
perovskite Sri - x Ba x Mn0 3 . Phys. Rev. Lett. 107, 137601 (2011). 

10. Giovannetti, G. et at Microscopic origin of large magnetoelectric coupling in 
Sr 1/2 Ba 1/2 Mn0 3 . Phys. Rev. Lett. 109, 107601 (2012). 

11. Amador, J. et at. Evidence of low-dimensional antiferromagnetic ordering and 
crystal structure in the R 2 BaNi0 5 (R = Y,Er) oxides. Phys. Rev. B 42, 7918 (1990). 

12. Darriet, J. & Regnault, L. P. The compound Y2BaNi05: A new example of a 
Haldane gap in a S — 1 magnetic chain. Solid State Commun. 86, 409-412 (1993). 

13. Alonso, J. A. et at Neutron diffraction study of the magnetic structure of 
Er 2 BaNiO s . Solid State Commun. 76, 467-474 (1990). 

14. Castro, M., Burriel, R., Salinas-Sanchez, A. and Saez-Puche, R. Heat-capacity 
measurement of R 2 BaNi0 5 (R = Y, Ho, Er) oxides. /. Magn. Magn. Mater. 
104-107,619-620 (1992). 

15. Garcia-Matres, E., Martinez, J. L. & Rodriguez- Carvajal, J. Neutron diffraction 
study of the magnetic ordering in the series R 2 BaNi0 5 (R — Rare Earth). Eur. 
Phys. J. B 24, 59-70 (2001). 

16. Garcia-Matres, E., Garcia-Munoz, J. L., Martinez, J. L. & Rodriguez- Carvajal, J. 
Magnetic susceptibility and field-induced transitions in R 2 BaNi0 5 compounds (R 
= Tm, Er, Ho, Dy, Tb, Gd, Sm, Nd, Pr). /. Magn. Magn. Mater. 149, 363-372 
(1995). 

17. Nenert, G. & Palstra, T. T. M. Magnetic and magnetoelectric properties of 
Ho 2 BaNi0 5 . Phys. Rev. B 76, 024415 (2007). 

18. Zheludov, A., Tranquada, J. M., Vogt, T. & Buttrey, D. J. Magnetic gap excitations 
in a one- dimensional mixed spin antiferromagnet. Nd 2 BaNi0 5 . Phys. Rev. B 54, 
7210 (1996). 

19. Yokoo, T., Zheludev, A., Nakamura, M. & Akimitsu, J. Ni-chain gap excitations in 
(Nd x Y!_ x ) 2 BaNi0 5 : One-dimensional to three-dimensional crossover. Phys. Rev. 
B 55, 11516 (1997). 

20. Maslov, S. & Zheludev, A. Coexistence of Haldane-gap excitations and long-range 
antiferromagnetic order in mixed-spin nickelates R 2 BaNi0 5 . Phys. Rev. B 57, 68 
(1998). 

21. Chen, J. W., Rao, G. N. & Li, K. W. Giant dielectric response of Haldane gap 
compound Y 2 BaNiO s . /. Appl. Phys. Ill, 064111 (2012). 

22. Singh, K. et al. Magnetoelectric coupling in the Haldane spin-chain system 
Dy 2 BaNi0 5 . Phys. Rev. B 88, 094438 (2013). 



23. Chepurko, G. G. et at Magnetic and spectral studies of Er 2 BaNi0 5 and other rare 
earth nickelates. Phys. Lett. A 157, 81 (1991). 

24. Basu, T., Iyer, K. K, Singh, K. & Sampathkumaran, E. V. Novel dielectric 
anomalies due to spin- chains above and below Neel temperature in Ca 3 Co 2 0 6 . 
Sci. Rep. 3,3104 (2013). 

25. Klimin, S. A., Kuzmenko, A. B., Popova, M. N., Malkin, B. Z. & Telegina, I. V. 
Infrared study of lattice and magnetic dynamics in a spin- chain compound 
Gd 2 BaNiO s . Phys. Rev. B 82, 174425 (2010). 

26. Sulewski, P. E. & Cheong, S-W. Raman-active phonons of the spin-1 chain 
antiferromagnets R 2 BaNi0 5 (R = Gd and Y). Phys. Rev. B 50, 551 (1994). 

27. Sulewski, P. E. & Cheong, S-W. Raman study of one-dimensional spin fluctuations 
in the spin-1 Heisenberg- chain compound Y 2 BaNi0 5 : Contrast with simple 
expectations. Phys. Rev. B 51, 3021 (1995). 

28. Fasolino, A., Los, J. H. & Katsnelson, M. I. Intrinsic ripples in graphene. Nat 
Mater. 6, 858-861 (2007). 

29. Perdew, J. P., Burke, K. & Ernzerhof, M. Generalized Gradient Approximation 
Made Simple. Phys. Rev. Lett. 77, 3865 (1996). 

30. Wu, H. et at Ising Magnetism and Ferrolectricity in Ca 3 CoMn0 6 - Phys. Rev. Lett. 
102, 026404 (2009). 

31. Choi, Y. J. et ah Ferroelectricity in an Ising Chain Magnet. Phys. Rev. Lett. 100, 
047601 (2008). 

32. Jahn, A. & Teller, E. Stability of polyatomic molecules in degenerate electronic 
states. I. orbital degeneracy. Proc. R. Soc. London Ser. A. 161, 220-235 (1937). 

33. Bersuker, I. B. The Jahn-Teller Effect and Vibronic Interactions in Modern 
Chemistry. (Plenum Press, New York, 1984). 

34. Albright, T. A., Burdett, J. K. & Whangbo, M.-H. Orbital Interactions in Chemistry 
(Wiley, New York, 1985). 

35. Pearson, R. G. The second-order Jahn-Teller effect. /. Mol. Structure: 
THEOCHEM 103, 25-34 (1983). 

36. Hippel, A. V., Breckenridge, R. G., Chesley, F. G. & Tisza, L. High dielectric 
constant Ceramics. Ind. Eng. Chem. 38, 1097-1109 (1946). 

37. Cohen, R. E. Origin of ferroelectricity in perovskite oxides. Nature 358, 136-138 
(1992). 

38. Maiti, K, Mahadevan, P. & Sarma, D. D. Evolution of electronic structure with 
dimensionality in divalent nickelates. Phys. Rev B 59, 12457 (1999). 

39. Kresse, G. & Hafner, J. Ab initio molecular dynamics for liquid metals. Phys. Rev. B 
47, 558-561 (1993). 

40. Kresse, G. & Furthmuller, J. Efficiency of ab-initio total energy calculations for 
metals and semiconductors using a plane-wave basis set. Comput. Mat. Sci. 6, 
15-50 (1996). 

41. Rodriguez -Carvajal, J. Recent advances in magnetic structure determination by 
neutron powder diffraction. Physica B: Condensed Matter 192, 55-69 (1993). 



Acknowledgments 

The authors thank Mr. S. Chowki for participation in sample preparation and Mr. Kartik K 
Iyer for his help in all these studies. 

Author contributions 

T.B. has carried out the magnetization, heat capacity, complex dielectric, resistivity and 
polarization experiments and analysed the data. B.G. and K.S. helped in polarization 
experiment and analysis. K.S. and N.P.L. carried out "/"-dependent XRD studies and S.B. 
carried out analysis of these XRD data. V.V.R.K. and P.M. has done the theoretical studies. 
S. Gohil and S. Ghosh carried out the Raman studies and analysed the data. N.M. and T.B. 
prepared and characterized the sample. E.V.S. has supervised the whole project, written the 
draft and finalized in consultation with other authors, in particular after discussions with 
P.M., S.G. (2), T.B. and K.S. 

Additional information 

Supplementary information accompanies this paper at http://www.nature.com/ 
scientificreports 

Competing financial interests: The authors declare no competing financial interests. 

How to cite this article: Basu, T. et al Displacive-type ferroelectricity from magnetic 
correlations within spin-chain. Sci. Rep. 4, 5636; DOI:10.1038/srep05636 (2014). 



| This work is licensed under a Creative Commons Attribution 4.0 International 
License. The images or other third party material in this article are included in the 
article's Creative Commons license, unless indicated otherwise in the credit line; if 
the material is not included under the Creative Commons license, users will need 
to obtain permission from the license holder in order to reproduce the material. To 
view a copy of this license, visit http://creativecommons.Org/licenses/by/4.0/ 



SCIENTIFICREPORTS | 4:5636 | DOI: 10.1038/srep05636 



6 



